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Because currently the most popular method of calculating plasma self- 
consistent fields is an incomplete method which, strictly speaking, is not suit- 
able to scientific investigation, we develop this method into be a complete 
reliable basic tool for scientific investigation. 
PACS: 52.65.-y, 52.35.-g. 
Plasma physics is a physical branch about many charged particles interacting through 
their self-consistent fields. In its earlier developing stage (about 1940s" 1960s), many the- 
oretical methods which are successful in other elder physical branch such as neutral gas 
physics and fluid mechanics were transplanted into this younger branch and rapidly built 
up the basis of this new branch. For example, in almost all basic textbooks in plasma 
physics [1-6], plasma wave is studied by not only non-relativistic fluid dynamics but also 
various analytical ansatzs on microscopic Vlasov- Maxwell equations. However, almost no 
one doubts whether these transplanted methods are appropriate for plasmas where numer- 
ous charged particles are correlated through their self-consistent fields. More important, 
in above-mentioned transplanted methods the plasma self-consistent fields is never strictly 
calculated but is indeed treated by various (obvious and hidden) approximations. 

Following example illustrates clearly a typical inconsistency of above-mentioned analyt- 
ical ansatzs on microscopic Vlasov-Maxwell equations. In many basic textbooks, people 
make Fourier analysis: fi = J^ fc fk exp(i6); E\ = Ek exp(i#);and 9 = kr — out on Vlasov 
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equation+Maxwell equations (in B = case) 

= d t f 1 +vVf 1 -E 1 -d p f ; (l.a) 

V • E 1 = J Ud 3 p (l.b) 

dtEx = [ vhd'v. (l.c) 



Especially, when Ei is of a mono chromic wave form E\ = Cons * exp [i (kr — ut)], this will 
lead to well-known Landau damping which is an important conception for plasma wave [1-6]. 
On the other hand, if E\ = Cons * exp [i (kr — ut)}, Eq.(l.b,c) will lead to 

vhd'v = ^ J hd'v 

and the fluid velocity meeting u = ^jj^ = f * JJ[^{^jJ^v - However, after deriving fluid 
momentum equation from Vlasov equation according to standard procedure, one could find 
following relation 



dt 



u J fid 3 v 



J fid'v + convective. term = E\ * J fod'v. (2) 



k ff 1 d*v + ff d*v_ 

This suggests that if the phase velocity f is a constant, the shape of the wave, or the shape 
of n = J fid 3 v, have to meet two equations, the continuity equation and Eq.(2) at d t ^ = 
case, and an assumed shape might not meet these two equations, which are both able to be 
expressed completely in term of n = J f\d 3 v. Therefore, a reliable investigation on a plasma 
wave should be based on strictly calculating self-consistent fields rather than some empirical 
ansatzs, which might be inconsistent to studied model equations. 

Although strictly calculating self-consistent fields represents a correct direction in plasma 
physics, the detailed technical road to achieve this goal causes this correct goal being greatly 
discounted. The chosen detailed technical road is the well-known particle simulation scheme, 
in which particles' information and fields' information are alternatively updated through nu- 
merous Newton equations and Maxwell equations [7-9] . This inevitably encounters a realistic 
question: the number of realistic particles is an astronomical figure and hence corresponds 
to too huge data mount. It is impossible to update so huge mount of data even one time. 



Therefore, this inevitably causes a merging approximation, which means TV realistic particles 
being merged into a so-called macroparticle and hence cut down significantly corresponding 
data mount, being introduced into the particle simulation scheme. Unfortunately, although 
this merging approximation makes updating being feasible, it will also cause a hidden dis- 
aster which refers to calculated E and B are functions of the merging ratio R me r g e = N : 1. 
People of course wish that with TV decreasing, the dynamics of calculated macroparticles is 
more and more approaching to the dynamics of real particles. Therefore, according to this 
viewpoint, it is reasonable to take the dynamics of calculated macroparticles of sufficiently 
small size as the dynamics of realistic particles. Unfortunately, this viewpoint is not true. 
It is almost impossible to give a mathematically strict proof on the uniform convergence 
of E (r,t; Rmerge — N : 1) to E (r, £; Rmerge = 1 '■ 1) with R m erge decreasing. Let us see an 
equation set including 2N + 4 equations describing N realistic particles (or macroparticles) 



d t E(R, t) = VxB(R,t) + J2 d ^ (t) 6 {n (t) - R) ; (3.1) 

i 

d t B(R, t) = -V x E(R, t); (3.2) 
V-E(R J t) = J2Hr i (t)-R); (3.3) 

i 

V-B(R,t) = 0; (3.4) 

dt , dtri(t) - = E(r t (t) , t) + dtn (t) x B(r t (t) , t) (3.2i+4) 

V 1 - Wi {t)f 

Vi = d t Ti (3.2i+5) 



If we merge two real particles into a macroparticle (or two macroparticles into a macro- 
macroparticle), we indeed deal with another equation set also including 2N + 4 equations 

d t E(R, t) = WxB(R,t) + ^2 dt r i (*) S (n (*) - R ) I 

i 

d t B(R, t) = -V x E(R, t); (4.2) 
V-E(R,t) = Y,Hr t (t)-R); (4.3) 



V-B(R,t) = 0; (4.4) 



d t r t it) 



V 1 - Wi it)f 

v-ik = d t r 2k 
d t r 2 k+i (t) , 



E{r t {t),t) + dtri (t) x B{n (t),t) 



d t r 2 k (t) 



\]l-[d t r 2k+1 (t)f \Jl-[d t r 2 k (t)f 
v 2 k+i = d t r 2k+1 



(4.4k+4) 

(4.4k+5) 
(4.4k+6) 

(4.4k+7) 



Here, if i is even or = 2k, Eq.(3.2i+4) and Eq.(3.2i+5) are just Eq.(4.4k+4) and Eq.(4.4k+5). 
But if % is odd or = 2k + 1, Eq.(3.2i+4) and Eq.(3.2i+5) will be replaced by Eq.(4.4k+6) 
and Eq.(4.4k+7). Obviously, the difference between these two equation sets is significant 
(because nearly half of total equations are approximated) and their respective solutions also 
have marked difference. More important, no matter how large N is, the difference between 
above two equation sets, and that between their solutions, is still significant. Therefore, it 
is impractical to expect following relation being valid 

d N I sup \E (r, t; R merge = N : 1) - E (r, t; R merge = 1 : 1) | J > (5) 

This relation, which implies that the smaller N is, the closer E (r,t;R merge = N : 1) is to 
E (r, t; Rmerge = 1 : 1), is just the theoretical basis of above-mentioned viewpoint. 

However, it is easy to obtain E (r,t; R merge = 1 : 1) if we note a fact: Any solution of 
following equation set of 2N + 5 members 







dt . "M*)»*) = _ E (r, it) ,t)-u (r, (t) , t) x B{ n (t) , t) 



y/l-[u(n (t),t)f 



n (6.o) 



d t E(R, t) = VxB(R,t) + Y^ d t r % (t) 8 (r ( (t) - R) ; (6.1) 

i 

d t B(R, t) = -Vx E(R, t); (6.2) 

V-E(R,t) = J2Hri(t)-R); (6.3) 



V • B(R, t) = 0; (6.4) 



d 



d t r l (t) u (n (t) , t) 



n = [d t r t (t) - u (r 4 (t) , t)} x B(r t (t) , t) (6.2i+4) 



l - [dtn (t)} 2 1 Ji-[u (n (t),t)} 2 



Vi = d t Ti 



(6.2i+5) 



... , 



where u (R, t) = £ 



i&n(t)=R 



dtr, (t) I £ 



i&i(t)=R 



1, must also be a solution of Eqs.(3). In a 



mathematical language, Eqs.(3) and Eqs(6) have their respective solution sets: {solutions 
of Eqs.(6)} and {solutions of Eqs.(3)}, and there strictly exists a relation between these two 
sets: {solutions of Eqs.(6)} C {solutions of Eqs.(3)}. On the other hand, one can easily find, 
by subtracting any two Newton equations in Eqs.(3), following relation for any two particles 



If two particles are at a same space position and the velocity of one particle is just equal 
to the fluid velocity at this position, Eq.(7) will automatically return to Eq.(6.2i+4) or 
Eq.(6.2j+4) at II = case. This fact implies that any solution of Eqs.(3) also meets Eqs.(6) 
at II = case, i.e., {solutions of Eqs.(3)} C {solutions of Eqs.(6) at II = case}. Because 
of the relation {solutions of Eqs.(6) at II = case} C {solutions of Eqs.(6)} C {solutions of 
Eqs.(3)}, thus, we have {solutions of Eqs.(3)} = {solutions of Eqs.(6) at II = case}. 

Eqs.(6) at II = case, a closed fluid equation set of u, E and B, implies a fast 
and exact method of calculating E (r,t; R merge = 1 : 1). Indeed, if people had noted fol- 
lowing strict relation several decades ago, the chosen detailed technical road will be free 
from above-mentioned disaster. Let us see well-known equation for fluid momentum 




E{r, (t) , t) + ckn (t) x B{r, (t) , t) - E{ rj (t) , t) + d t r 3 (t) x B(r 3 (t) , t). 



(7) 



Pfi = Jpfd 3 p/ J fd 3 p [1-6] 



dtPfi + Ufi ■ Vp/z = E + Ufi x B + thermal pressure/density 



(A.l) 



where Ufi is the fluid velocity Ufi — J vfd 3 p/ J fd 3 p. It is a pity for the community of 



5 



plasma physics during several decades to fail of finding a very simple short-cut from above 

equation. Because the velocity v is a nonlinear function of the momentum p (i.e., v = J- — ) 

V 1+ p 2 

and vice versa, we should note that the statistic average value J pfd 3 p/ J fd 3 p (i.e. fluid 
momentum) is usually not equal to the momentum corresponded by the statistic average 
value J vfd 3 p/ J fd 3 p (or fluid velocity), i.e., pfi ^ p (ufi)(where p(ufi) refers to the value of 
function p (variable) at variable = Ufi), if the distribution / is not a Dirac function of p (i.e., 
/ has a thermal spread over p-space). Only at zero temperature case, there is p^ = p(ufi). 
(Strictly speaking, if / is a symmetric function of p, there will be pfi = p(ilfi) = 0, Ufi = 
and thermal pressure^ 0. But this special case corresponds to E = and E + Ufi x B = 0. A 
non-zero thermal pressure will drive pfi differing from according to Eq.(A.l). Once pfi ^ 0, 
there will be pfi ^ p(ufi) because / has an asymmetric thermal spread over p-space). 

Any distribution function / has two independent characteristic parameters: the variance 
and the mean. Here, the mean of / is represented by Uf\. As above discussed, Pfi is thus 
a binary function :pfi = pfi (uf i, variance) . Thus, we could express pfi as a series of vari: 
Pfi = J2 i>0 dl ari pfi\ vari=0 (vari) 1 . Moreover, a term in Eq.(A.l), thermal pressure/density, 
could also be expressed as a similar series of vari but this series does not contain (uari)°-term 
because when vari = 0, / is a Dirac function and hence thermal pressure is equal to zero. 
Namely, thermal pressure/density= Y2i>i c * ( var i)' ' ■ Note that Eq.(A.l) is valid for any value 
of vari. Then by substituting these two series into Eq.(A.l) and comparing coefficients of dif- 
ferent order (vari) 1 -terms, we could find an equation set of infinite members and every mem- 
ber is for a definite order (uari) J -term (because the fact that Eq.(A.l) is valid for any value 
of vari requires these coefficients of different (vari )*-terms being zero separately). For zero- 
order term ( (vari )°-term), we have = d t \pfi (u^, 0)] + Ufi ■ V r [p/; (u^, 0)] — [E + x B]. 
Indeed, all equations for non-zero-order terms ((vari) 1 ^ -terms) could be merged into an 
equation if we substracting Eq.(A.l) and the equation for zero-order term ((vari) -term), 
and this equation reads = d t [pfi — Pfi (m/z, 0)] + Ufi ■ V r [Pfi — Pfi (ufi, 0)]+thermal pres- 
sure/density. Note that pfi (ufi,0) = p(ufi) = , Ufl 2 . Namely, because Eq.(A.l) is of a 
general form = function(ufi,vari), the conditions for such a general form being valid at 



arbitrary vari-valne are: = function(ufi,0) and d var i function = 0. For Eq.(A.l), these 
two conditions could be expressed by two equations 

= dt\pfi - Pfi (uf h 0)] + Ufi ■ V r [Pfi - Pfi ( u fh 0)] + thermal pressure/density; (A. 2. a) 
= d t \p fl {u fl , 0)] + ufi ■ V r [p fl (u fh 0)] - [E + u fl x B] , (A.2.b) 

and the former equation depends on the temperature (or vari) whereas the latter is inde- 
pendent of the temperature. Obviously, the solution of Eqs.(A.2) must be the solution of 
Eq.(A.l). Thus, Ufi, E and B form a closed description. 

Obviously, the closed equation set (Eq.(A.2.b)+4 Meqs) has been misjudged as only valid 
at zero-temperature case for a long time, and therefore is often judged as an "approximation 
theory". People are easy to believe that the thermal pressure will inevitably affect pfi and 
"further affect self-consistent fields (through Meqs)". Actually, a scrupulous reader will not 
have such a misunderstanding because he could notice that it is Ufi, rather than pfi, that 
appears in Meqs. Moreover, Eq.(A.2.a) will lead to infinite correlated quantities which are 
all dependent on the variance of /. Like the treatment in many standard textbooks [1- 
6], introducing an assumed thermodynamics state equation could cut the chain of infinite 
correlated quantities and form a closed equation set. But this will again depreciate greatly 
the validity of fluid theory. Actually, fluid theory could become a nearly perfect basic 
tool in plasma physics according to following way: First, figuring out E and B of realistic 
particles and Ufi through Eq.(A.2.b)+4 Meqs. Second, if one is interested in more detailed 
particles' information (such as particles trajectories and velocities, or so-called phase-space 
snapshots), he could solve Vlasov equation under known E and B, or make a "test-particle 
calculation" on macroparticles, and then extract relevant information from solved / or solved 
trajectories and velocities of macroparticles. More important, it is no need to introduce an 
assumed thermodynamics state equation. Eq.(A.2.a) itself actually implies a complicated 
thermodynamics state equation. 

A more direct way of obtaining above closed equation set of Ufi, E and B could 
be start from Vlasov equation (VE). For any distribution /, we could expressed 
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as / = nS (v - u f i) + a S (v - u f i) + Y,i a z>i( v ~ u fiY (where n = J fd 3 p, u fl = 
J vfd 3 p/ J fd 3 p, cii are independent of v, a® depends on all coefficients a^i through a 
relation J a 6 (v — Ufi) + X^i>i a « ( v ~ u fiY d 3 p = 0, i.e. ao is a function of all a^i, 
a o — Oo (ai, di, ...)). Substituting this expression into VE and comparing the coeffi- 
cients of (v — UfiY-tevm, we could find that there exists following equation for / mono = 
nb (v — Ufi) + ao6 (v — Ufi) (because of the fact that VE is valid for any distribution / which 
is characterized by its own coefficient set {a«}) 



dtfrnono + Ufl ■ V ' frnono ~ [E + Ufi X B] ■ d p f v 



0. 



(A.3) 



which could directly lead to Eq.(A.2.b) according to standard procedure. 

Moreover, there is an easier way of obtaining this closed equation set of Ufi, E and B 
from the starting equations of particle simulation. Note that a relativistic Newton equation 

dtn it) 



= d t 



-E(n (t),t)-dtri (t) x Bin {t),t) 



u (n (t) , t) 



dtTj (t) 

u (n (t) , t) 



[d^ (t) - u fa (t),t)} x B(r t (t) ,t) 



d 



- E fa (t) ,t)-u fa (t) , t) x B(n (t) , t) 



(B.l) 



'l-[u(r, (i),i)] 2 

is valid for arbitrary value of dtVi (i), or arbitrary value of A = c^r, (t) — u fa fa , t). Because 
Eq.(B.l) is of a general form = function(u fa (t) , t) , A = aVi (i) — m (rj (t) , t)), the con- 
ditions for such a general form being valid at arbitrary A-value are: = function(u, 0) and 
d&f unction = 0. For Eq.(B.l), these two conditions could be expressed by two equations, 











d t r l (t) 



d f 



u (n fa , t) 



dt 



u (^ fa , t) 



[d^fa-ufa (t),t)] x 5fa fa fa 



(B 



- £ (rj fa , t) - u (rj fa , t) x 5fa fa , t) 



:b 



'l-[u(r, fa fa] 2 

where w (i?, t) represents the average value of the velocities of macroparticles whose positions 
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are at R when the time is t, i.e., u(R,t) = J2i^ n (t)=R ^t r i W / J2i£ n (t)=R 1- Obviously, 
Eq.(B.2.b) is just equal to Eq.(A.2.b). 

This also implies that even in the particles simulation scheme, it is no need to alterna- 
tively update particles' data and fields' data because fields' data could be strictly updated 
through Eqs. (6.0-5) at II = case. Here, alternative updating is the basic reason for the 
particles simulation scheme being entangled with above-mentioned disaster. This elder road 
characterized by alternatively updating, is indeed an approximation on the starting model 
equations Eqs. (3). If this disaster is not overcome, strictly speaking, particle simulation 
scheme is an incomplete method not suitable to scientific investigation. Here, this disaster 
could be easily overcome if another detailed road, which does not involve any approximation 
on Eqs. (3), is adopted, i.e., directly solving E and B from Eqs. (6.0-5) at n = case, rather 
than alternatively updating fields and particles. In short, no matter which one of the par- 
ticles simulation scheme, the relativistic fluid theory, and the microscopic Vlasov-Maxwell 
theory is chosen by people as the starting model of investigating plasma physics, the self- 
consistent fields, E and B, obey a fixed fluid equation set, i.e., Eqs(6.0-5) at n = case or 
Eq.(A.2.b)+4 Meqs. Indeed, these different basic methods are equivalent if they are in their 
respective strict forms. There is no reason to think that any method is better than others. 

A more concise and straightforward presentation of above text after the paragraph around 
Eq.(7) is given as below: 
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Eqs. (6.0-4) at II = case, a closed fluid equation set of u, E and B, implies a fast and 
exact method of calculating E(r,t; R mer ge = 1:1). 

Indeed, if we notice the universality of some physical laws, we could obtain this closed 
equation set of u, E and B from different basic theoretical methods in plasma physics. 
At present, there are mainly three basic theoretical methods: fluid theory, Vlasov-Maxwell 
(VM) theory and above-mentioned particle simulation (or well-known PIC method). In every 
basic theoretical method, except 4 Maxwell equations (Meqs), other equations correspond 
to a basic physical law. 

l.In fluid theory, this basic physical law is represented by fluid momentum equation [1-6] 

= d t pfi + Ufi ■ Vpfi — E + Ufi x B — thermal pressure/density 
= [dt(j>fi — p(ufi)) + Ufi ■ V (pfi — p{ufi)) — thermal pressure/density] 

+ [dtp (u n ) + u n ■ Vp (u f i) -E + u fl xB}; (A.l) 

where Ufi is the fluid velocity Ufi = J vfd 3 p/ J fd 3 p and Pfi = J pfd 3 p/ J fd 3 p is the 
fluid momentum. Because the velocity v is a nonlinear function of the momentum p (i.e., 



v 



and vice versa, we should note that the statistic average value J pfd 3 p/ J fd 3 p 



(i.e. fluid momentum) is usually not equal to the momentum corresponded by the statistic 
average value J vfd 3 p/ J fd 3 p (or fluid velocity), i.e., pfi ^ p (ufi)(where p(ufi) refers to 
the value of function p (variable) at variable = Ufi), if the distribution / is not a Dirac 
function of p (i.e., / has a thermal spread over p-space). Only at zero temperature case, 
there is pfi = p(v,fi). (Strictly speaking, if / is a symmetric function of p, there will be 
Pfi = P( u fi) = 0; u fi = an d thermal pressure^ 0. But this special case corresponds to 
E = and E + UfixB = 0. A non-zero thermal pressure will drive p^ differing from 
according to Eq.(A.l). Once pfi ^ 0, there will be pfi ^ p(ufi) because / has an asymmetric 
thermal spread over p-space). 

2. In VM theory, this basic physical law is represented by Vlasov equation [1-6] 

= d t f + vVf-[E + vxB]-d p f. 
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[dt(f ~ fmono) + V ■ V (/ - fmono) — [E + V X B] • d p (f — fmono)} 

+ (V -Ufl) ■ V f mono - (V -Ufl) ■ d p f mon o 

+ [dtfmono + Ufl ■ V fmono — [E + Ufi X B] ■ dpfmono] 



(B.l) 



Any distribution function / has two independent characteristic parameters: the variance 
and the mean. Here, the mean of / is represented by Ufi = J vfd 3 p/ J fd 3 p). For any 
distribution /, we could expressed as / = n5 (v — Ufi) + a 6 (v — Ufi) + Yli>i a i ( v ~ u fiY 
(where n = J fd 3 p, Ufi — f vfd 3 pj J fd 3 p, are independent of v, a depends on all 
coefficients through a relation J a^b (v — Ufi) + X^>i a i ( v ~ u fiY d 3 p = 0, i.e. clq is a 
function of all a^i, a$ = a (a\, a iy ...)). The variance, vari, depends on the values of all 

3. In particle simulation, this basic physical law is represented by relativistic Newton 
equation of any particle [7-9] 

d t ri it) 



= d t 



d 



E(n (t),t)-dtri (t) x B( ri (t),t) 

u (n (t) , t) 



d 



1 - [dtr, (t)} 2 
dtn (t) 

' ^i-idknit)) 2 1 ^i-iu( ri (t),t)] 

U (Ti (t) , t) 



dt 



[d t n (t)-u(n (t),t)] x B{r % {t),t) 



- E (ri (t) ,t)-u (r { (t) , t) x B(n (t) , t) 



(C.l) 



'l-[u(r, (t) } t)Y 

where u (R, t) represents the average value of the velocities of macroparticles whose positions 
are at R when the time is t, i.e., u (R, t) = Y^i &i {t)=R (*) / J2ien(t)=R 1- 

The common of Eq.(A.l), Eq.(B.l) and Eq.(C.l) is that they have a binary- function type 
general form 



= function(varl,var2) 



(8) 



where varl and var2 are independent variables. 1. For fluid theory, because both 
Pfi — p(ufi)and thermal pressure/density depend on vari (i.e., the variance of /), there are 
varl = Ufi and var2 = vari (i.e., the variance of /). Note that thermal pressure/density = 
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^2i>i c i { var "iY because when vari = 0, / is a Dirac function and hence thermal pres- 
sure is equal to zero and hence this series does not contain (varif '-term. Likewise, pfi 
is thus a binary function :pfi = pfi (ufi,vari) and could be expressed as a series of vari: 

Pfi - Y.i>o dlariPfi\vari=o (vari) 1 and pfi (ufi, 0)=p (ufi) = -7=5= ■ 2. For VM theory, there 

V u fi 

are vari = f m0 no and var2 = v—Ufi. 3. For particle simulation, there are vari = u (r, (i) , t) 
and var2 = d t ri (t) — u (rj it) , t). 

The universality of these basic physical laws requires that they are valid for arbitrary 
value of var2. For example, one cannot expect that a particle does not meet relativistic 
Newton equation when it corresponds to var2 = d t ri (t) — u (r, (t) , t) — (because following 
equation is absurd) 

dt 1 - £7 (r t (t) ,t)-u {r t (t) , t) x ^(r, (t) , t) ? 0. 

Likewise, following equation is also absurd 

[%> (w/i) + u/i • Vp - E + u n x B] 7^ 0. 

because it implies that the fluid momentum equation is invalid at zero-temperature var2 = 
vari = 0. 

The conditions for such a universality requirement are: 

= function(var 1,0); (91) 
= d var2 function. (9.2) 

Let us see detailed forms of these two conditions. 1. For Eq.(A.l), these two conditions 
could be expressed by two equations 

= d t [p f i - p f i (uf h 0)] + u f i ■ V r [Pfi ~ Pfi {ufi, 0)] + thermal pressure/density; (A. 2. a) 
= d t [p fl {u fl , 0)] + ufi ■ V r [p fl (u fh 0)] - [E + u fl x B] . (A.2.b) 

2. For Eq.(B.l), we have 
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= [8t(f - fmono) + V ■ V (/ - fmono) ~ [E + V X B] • 8 p (f - fmono)} 

+ (V - Ufl) ■ V f mon o - (V - Ufl) • d p f mono \ (B.2.a) 
= [3 t fmono +Ufr V fmono ~ [E + Ufl X B] ■ d p fmono] ■ (B.2.b) 



3. For Eq.(C.l), we have 



d* 



u (r l (t) , t) 







dt 



u (n (t) , t) 



[dkn {t)-u(n (t),t)} x B{n (t) ,t) 



(C.2 



- E {n (t) ,t)-u (n (t) , t) x B{n (t) , t) 



'i-[u( ri (t),t)] 2 

Obviously, Eq.(A.2.b) and Eq.(C.2.b) are same, and Eq.(B.2.b) could also be trans- 
formed into them according to standard procedure. Therefore, Eq.(A.2.b) or Eq.(B.2.b) 
or Eq.(C.2.b)+4 Meqs forms a closed equation set. 



(C.2 
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